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Abstract 



I present a short review of irradiation processes in close binary systems. 

1. Introduction 

Thermal- viscous instability is widely accepted as the origin of Dwarf Nova 
and Soft X-ray Transient (X-ray Nova) outbursts. In the 'standard' model, one 
assumes a constant mass-transfer rate, no effects of irradiation are taken into 
account and the accretion disc is supposed to extend down to the surface of the 
accreting body (or the last stable orbit in the case of accreting black holes and 
compact neutron stars). This version of the model, however, explains neither 
all the observed outburst types nor all the observed outburst properties. Gener- 
alizations of the 'standard' model which take into account some or all of these 
neglected, but obviously important, effects have been proposed and studied in 
various investigations. I will shortly discuss some recent results concerning irradi- 
ation (Hameury, Lasota & Hure 1997; Dubus et al. 1999; hereafter DLHC; King 
(1997); King & Ritter 1998). 

Since in the literature of the subject most articles use incorrect equations 
to describe irradiated accretion discs I begin by repeating the simple derivation 
of the basic equations presented in DLHC. 

2. A simple introduction to the vertical structure of irradiated discs 

In accretion discs the vertical energy conservation equation has the form: 
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Fig. 1. Vertical structure of an accretion disc around a M = 10M Q compact 
object at r = 3 x 10 10 cm. a « 0.1, M « 10 16 g s" 1 (i.e. T efr « 5700 K) 
and Ti rr = 0. Both temperatures are normalized by 10 4 K. Since Tj rr = 0, the 
surface temperature T(r s ) = T e g. The dashed line is the density in units of 
xlO -7 g cm -3 , the dotted line the ratio of the convective to the total fluxes. 
Note that for these parameters, the section of the disc lays on the lower stable 
branch. (From DLHC) 



where F is the vertical (in the z direction) radiative flux and Q v i S (R,z) is the 
viscous heating rate per unit volume. Eq. ([!]) states that an accretion disc is 
not in radiative equilibrium, contrary to a stellar atmosphere. Using the "a - 
viscosity" prescription (Shakura & Sunyaev 1973) v = (2/3)ac^/flK, where a 
is the viscosity parameter (< 1), f2x is the Keplerian angular frequency and 
c s = \jPjp is the sound speed, p the density, and P the pressure one can write 

Q vis (R,z) = (S/2)an K P (2) 

Viscous heating of this form has important implications for the structure of opti- 
cally thin layers of accretion discs and may lead to creation of coronae and winds 
(Shaviv & Wehrse 1986; 1991). Here, however, we are interested in the effects of 
irradiation on the inner structure of an optically thick disc, so our results should 
depend on the precise form of the viscous heating. We neglect, however, the pos- 
sible presence of an X-ray irradiation generated corona and wind, described by 
Idan & Shaviv (1996). When integrated over z, the rhs of Eq. ([]]) is equal to 
viscous dissipation per unit surface: 

F vis = -aQ K / Pdz, (3) 
2 Jo 
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Fig. 2. Same as Fig. 1. but for T irr = 12000 K. (From DLHC) 

which is close, but not exactly equal, to the surface heating term (9/8)z/£fi|- 
generally used in the literature. The difference between the two expressions may 
be important in numerical calculation (Hameury et al. 1998) but in the present 
context is of no importance. One can rewrite Eq. ([[]) as 

d,F F ■ 

where we introduced a new variable, the optical depth dr = —n-Rpdz, /% being 
the Rosseland mean opacity and tt = / +o ° K^pdz is the total optical depth. As 
shown in DLHC, putting /(r) = 1 is a good approximation. 

At the disc midplane, by symmetry, the flux must vanish: F(tt) = 0, 
whereas at the surface, (r = 0) 

F(0) = aT e 4 ff = F vis (5) 

Equation (|) states that the total flux at the surface is equal to the energy dis- 
sipated by viscosity (per unit time and unit surface). The solution of Eq. (|J) is 
thus 

F{t) = F vis (l - ^) (6) 

where r tot is the total optical depth. 

To obtain the temperature stratification one has to solve the transfer equa- 
tion. Here we use the diffusion approximation 
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appropriate for the optically thick discs we are dealing with. The integration of 
Eq. (0) is straightforward and gives : 

T\r)-T\0) = 3 -r(l-^jT^ (8) 

The upper (surface) boundary condition is: 

T\Q) = \T^ + n r (9) 

where Tj 4 r is the irradiation temperature, which depends on R, the albedo, the 
height at which the energy is deposited and on the shape of the disc (see Eq. |13|) . 
In Eq. (H) T(0) corresponds to the emergent flux and, as mentioned above, T e s 
corresponds to the total flux, hence the factor 1/2 in front of T e g. The temperature 
stratification is thus : 

+ n r . (io) 

For 7t 3> 1, the temperature at the disc midplane is 

Tc =T 4 (tt) = jjr tot T e 4 ff + 7] 4 r (11) 

It is clear, therefore, that for the disc inner structure to be dominated by irradi- 
ation and the disc to be isothermal one must have 

F rxT 4 

£JE = » Fvis (12) 

r tot 7~tot 

and not just F irr 3> F v - m as is usually assumed. The difference between the two 
criteria is important in low-mass X-ra binaries since, for parameters of interest, 
rtot ~ 10 2 — 10 3 in the outer disc regions. 

The effect of disc irradiation is illustrated on Figs. 1 & 2 (DLHC). Fig. 1 
shows the vertical structure of a ring which is part of an unilluminated accretion 
disc. This ring is on the lower, cool branch of the S-curve. The energy transport 
is dominated by convection. The surface temperature is equal to the effective 
temperature given by viscous dissipation (see Eq. |5|). The vertical structure of 
an irradiated disc is shown on Fig. 2. Although the irradiating flux is 20 times 
larger than the viscous flux the disc is not isothermal. Since in the non-irradiated 
disc T c « 14500 K and the surface temperature is T s « 5700 the optical depth is 
~ 100 and, obviously an irradiation temperature higher than 14500 K is required 
to make the disc isothermal. Note that irradiation suppressed convection and the 
disc now is purely radiative. 
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Fig. 3. £ - T surface , E - T c , E - T eff and E - M curves for r = 3 • 10 10 cm, 
M = 1OM , a = 0.1, and T irr = [0, 3, 6, 9, 12] x 10 3 K. (From DLHC) 



3. Can the outer disc 'see' a point source located at the midplane ? 

The answer is 'no' (see e.g. DLHC). The reason is that contrary to the 
frequently made assumption a 'standard' accretion disc is convex rather than 
concave. Images showing flaring discs have nothing to do with reality, at least 
with reality described by the standard model of planar discs. 

For a point source, the irradiation temperature can be written as 

4 = V Mc*(l - e) H- ai ( d\nH ix 
irr AnaR 2 R { d\nR 

where n is the efficiency of converting accretion power into X-rays, e is the X-ray 
albedo and ifi rr is the local height at which irradiation energy is deposited, or the 
height of the disc "as seen by X-rays". We use here Hi rr and not H, the local 
pressure scale-height, as is usually written in the literature because, in general, 

Hi TT ^ H. 

Eq. (|13D is usually used (at least in the recent, abundant, publications on 
the subject) with 'typical' values of e > 0.9, H/R = 0.2 (no difference is seen 
between H iTT and if; the fact that there is no reason for the photosphere to be 
at one (isothermal) scale- height seems to be largely ignored). These values are 
supposed to be given by 'observations'. However, when one reads the articles 
quoted to support this assertion (Vrtilek et al. 1990; de Jong et al. 1996) one 
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Fig. 4. Radial profiles of the midplane temperature, surface density, and photo- 
spheric height to radius ratio for an un-irradiated (continuous line) and irra- 
diated disc (dashed line) around a 10 M compact object. In the temperature 
diagram the continuous line represents the effective temperature. The accre- 
tion rate is M = 10 18 g s -1 , T irr is taken from Eq. (15) with C = 5 x 1(T 4 . 
Regions beyond the radius at which a break in the T c or X curve appears are 
unstable. (From DLHC) 



finds nothing of the kind there. These papers assume that an irradiated disc is 
isothermal, de Jong et al. (1996) fit light-curves with the isothermal model of 
Vrtilek et al. (1990). Moreover, since de Jong et al. (1996) model lightcurves 
of neutron-star binaries the H/R = 0.2 cannot be applied to black-hole binaries 
(especially if H were the pressure scale- height): since black- holes are more massive 
than neutron stars, in their case H/R should be smaller at a given radius (let me 
add for the benefit of some readers: this is because gravity is then stronger), as 
seen in Fig. 6 & 7 of DLHC. In any case, for H/R = 0.2 the vertical hydrostatic 
equilibrium equation would imply high temperatures at the disc's outer rim: 

\M Q J V10 10 cm/ \RJ ' V ; 

which is clearly contradicted by observations. 

When one calculates self-consistent (in the sense that ifi rr in Eq. (13) is 
calculated and not assumed) models of irradiated discs one sees that the outer 
disc regions whose structure could be affected by irradiation are hidden in the 
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shadow of the convex disc. This 'self-screening' results from the same physical 
process that is at the origin of the dwarf-nova instability: a dramatic change of 
opacities due to hydrogen recombination. Therefore, although irradiation could 
stabilize the disc, the unstable disc regions (see Fig. 3), cannot be irradiated 
by a point source located at the disc mid-plane. Models invoking the stabilizing 
effects of irradiation (van Paradijs 1996; King & Ritter 1998) have therefore to be 
revised (DLHC). Since outer disc regions in low mass X-ray binaries are clearly 
irradiated, (van Paradijs & McClintock 1995) model revisions must concern the 
disc-irradiating source geometry. To represent a geometry allowing the disc to 
see the irradiating source DLHC assumed that 

4 Mc 2 

Tta = (15) 

and calculated irradiated disc structure. The results are shown in Fig. 4. The 
continous line represents both a non-irradiated disc and a disc irradiated according 
to Eq. (13) with a self-consistently calculated H irr : there is no difference between 
the two cases. 

4. Can irradiation by a hot white dwarf explain the UV-delay in dwarf 
novae? 

As shown by Hameury, Lasota & Dubus (1999) the answer to this question 
asked by King (1997) is 'no'. Irradiation of the disc by the hot white dwarf may, 
however, be important in a different context (see Sect. 6). 

5. Can the secondary in a dwarf-nova system be irradiated during 
outbursts ? 

The answer is 'yes' at least for SS Cyg (Hessman 1984) and WZ Sge (Smak 
1993). In these systems it was observed that the secondary's hemisphere facing 
the accreting white dwarf was, during the outburst, heated to 16 000 - 17000 K, 
which for WZ Sge implies that the irradiating flux is ~ 5000 times larger that the 
star's intrinsic flux. It is hard to believe that the companion does not increase its 
mass transfer rate in such conditions. WZ Sge is a very special system anyway 
(Lasota, Kuulkers & Charles 1999). 

6. Can irradiation of the disc and of the secondary determine outburst 
properties? 



Warner (1998) suggested that outburst properties of SU UMa stars could 
be explained by the effects of irradiation of both the accretion disc and the sec- 
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ondary. Preliminary results by Hameury & Lasota (1999, in preparation; see also 
Hameury these proceedings) seem to confirm that properties of SU UMa (and ER 
UMa) stars may be explained in this way. 
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